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Abstract 

A time dependent generalization of the Ginzburg -Landau Lagrangian is proposed. It contains two terms 
determining the time dependence and the four arbitrary scalar functions. Relevant equations, which coincide 
with equations following from the suitable Hamiltonian, are derived by a standard variational technique. These 
equations determine the energy conservation law and admit twofold time dependence which leads either to first or 
to second order time derivatives in Ginzburg - Landau equations. By introducing the gauge invariant potentials 
and choosing the gauge which differs slightly from the classical Lorentz one, the theory simplifies significantly. 
The results gained are discussed and compared to some earlier propositions. The presented approach, when 
reduced to a static one, is found to be in perfect agreement with that reported recently by the Kolacek group. 
This indicates indirectly, that the equation with the first order time derivative seems to be more justified. 

PACS numbers: 74.20.De, 74.50.+r 



1 Introduction 

The Ginzburg - Landau (GL) free energy and the derived then consequently GL equations have proven to be 
very fruitful in the theory of superconductivity. Different situations, a lot of extensions, generalizations and 
applications to the different type of superconductors can be found in a huge and commonly known literature 
(c.f. ^ , 1^). But, perhaps as a natural consequence of the original article the majority of papers and 
textbooks starts from a static version, when time is completely ignored. 

The aim of the present paper is to correct this deficiency and to incorporate time in the frame of GL 
formalism, so that the obtained relations can be useful in the analysis of dynamical problems in the theory of 
superconductivity. The static reduction should then coincide with the known relations, issuing from the free 
energy concept. 

There were different attempts to include time in the GL formalism. These attempts can be divided into two 
groups. In the papers of the first group, usually the time dependent formulation is postulated, starting a priori 

from the time-dependent Ginzburg - Landau (TDGL) equation ^hipt = {2m)~^ (fiV)^ + a — /S ip, as in Q 
or similar. Its applicability to superconductive problems is later discussed and investigated. We are aware of a 
rich literature devoted to the mathematical aspects of this equation and with a proper choice of constants with 
its affinity to the nonlinear Schrodinger equation, see e.g. Q. The second group, far less numerous, is formed 
by publications, in which the authors start by choosing some reasonable Lagrangian, as e.g.[^, j^]. There are 
also some papers, in which time is introduced from other premises e.g. 

We would like to point out that in the time dependence is involved also the problem of the scalar potential 
and charge density. Can the gauge be chosen such that scalar potential can vanish, and what is the role of a 
nonzero charge density then? In many textbooks considering static consequences of GL equations, the scalar 
potential is assumed to be zero without convincing arguments. Better insight into this problem can gained from 
references 0, and references therein listing earlier contributions. In this paper we propose a modification 
of the GL formalism in order to include the time dependence starting from a suitable Lagrangian. Since the 
authors are convinced that there may be different approaches to the problem, in the first part there is proposed 
a general form of a GL-like Lagrangian, in which time is included in two different manners. Moreover, there 
are included four arbitrary scalar functions, which give additional freedom and assure some universality of 
considerations. In principle, the whole approach is gauge invariant. Next, the relevant modified GL equations 
are derived. These, however, are rather complicated and illegible, when as independent quantities are chosen 
the standard order parameter and potentials. Nevertheless, the formalism can be considerably simplified by the 



introduction of so called gauge - invariant (g.i.) potentials. This concept is only partially new, since the first of 
these potentials is the commonly known gauge invariant phase. 

The usefulness of this approach for superconductivity is discussed in the next sections. We show that our 
treatment in the static case is in agreement with the static approach reported in references Q, 1^, More 
precisely, some particular static reduction of our general theory is in very good agreement with the references 
cited, with an accuracy to a single term, which exists in our theory, but probably was simplified in the cited 
papers. This result serves as an important hint concerning the form of TDGL. The discussed agreement indicates 
that the quantity, which has physical interpretation as charge density must be independent of scalar potential. 
A natural consequence is such reduction of TDGL that leads to the presence of the first order time derivative 
in GL equation , which however seems to be in contrast to the case considered in reference ^ . 

This last result automatically determines the final form of the TDGL-like Lagrangian and, as a consequence, 
the TDGL equations. 

2 Time dependent Ginzburg - Landau approach 

We propose the extended version of the Lagrangian density which includes a time dependence in the form 



GL 



(1) 



jC-gl = 



(At + V^)' - (V X A)' 



W ~V[{-iV - A)^jj] {iV~A)4> 



Ct = M i Udtil^ - ^pdti:] - 2^i^Lp + N[{idt~' {-idt-f^^p 



(2) 
(3) 

The canonical quantities here are: vector and scalar potentials A and (p, complex order parameters ip and 
ip, and the density of the charge carriers n. We assume that the functions W, V, M, N are algebraic real functions 
of two arguments: [ipl and n, i.e. W = W {q,n) with.g := lipl , etc. and play a role of weight functions which 
can be fixed for concrete problems. Later we will impose the symmetry condition that = iIj*. 

Choosing e.g. W — jV'l^ + 2 and dropping the term (Af + Vip)'^ , it is seen that £gl coincides with 
the standard, static version of GL free energy which can be found elsewhere, (e.g. |^ ). Some particular 
choices of W,V, M, N were reported in and ||], announced also in||ll[ and in the static version discussed 
in numerous books and papers e.g. [||, . Thus the Cgl describes the static part of the Lagrangian 
density and Ct, the supplement introducing a twofold time dependence in further field equations. We introduce 
canonical variable n describing the density of free charge carriers under infiuence of the lecture and paper p!o[ |. 
Of course all these quantities depend on space coordinates x, y, z and time t. 

The proposed Lagrangian density is gauge invariant with respect to the transformation 



{A: 



Vx, ^ + xt, -0 ^ f/'exp {ix) , -0 ^ V'Gxp (-ix) , 



' + X, n^n\, 



(4) 



which will be evident shortly below (under assumption = ip*). 

The standard variational procedure leads to the system of rather complicated equations 



Sa ■■ [At + Vip)t + V X A + V 
6^ : -Alp - V (At) - + N 



5n 



= 



Wn + Vn [(-iV - A) V] (jV - A) ^ 



'Nn [{idt - (p) tp] {-idt - (p)ip 



i {ipdtilj — ipdtifj — 2'ipipip 



= 



W.-V.[{tV + A)ij] (zV-A)^ 



+Nr [{tdt - ip) iP] m + 



V 



-Aip- 



i {%pdtip — TpBtTp^ — 2ip'ipip 
(V-V • A + 2A-VV') + A^V-] 



+N [df^p + i i^dtip + 2ipdt^p) - ip'^ip] + 2M [^Pip - idti^] 
- {Vip - iAiP) VV + {dtijj + iiptp) dtN - iqdtM = 



(5) 
(6) 

(7) 



(8) 
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Equations and can be rewritten as nonhomogeneous wave equations for potentials 



AA - Att = V \i - ipyipj + 2A?/;-0j + V (V • A + vJt) 

Aip - iptt = -^Mi^i) + N 2(p'iPi} + i (^V'V't - V'V't) - (V • A + (^t) 



(9) 
(10) 



The last terms in both equations vanish, if the Lorentz gauge V • A + (y9( = is assumed. Then the first 
terms on the right hand sides represent the current j and charge p densities, respectively. Thus 



i = V i (^tpVip - ipVifj + 2Aip-ip 



p = 2Mipip - N 



(11) 
(12) 



Assuming ijj — tp* and calculating the imaginary part of [^) exp {—iO}] one obtains immediately the 
continuity equation V • j + pt = with j and p given by ([ll| ) and (12), respectively. 

Thus in the presented formalism we are able to include time dependence and four different types of weight 
functions. The physical interpretation will be evident later on. Here, it is worthwhile to remark that a uniform 
and time independent charge density of the background or lattice can be introduced by a suitable choice of M 
function. 

As it was mentioned above, all derived relations are gauge invariant with respect to (^), although it is not 
seen directly. It can be done quite immediately defining new variables F and / which surely are gauge invariant 
by relations 



F:=V9-A, f:=^- 



(13) 



It is obvious that the first quantity F is the gauge invariant phase commonly introduced description different 
superconductive problems. The second one, / was to our knowledge never defined. Both quantities redefine the 
vector and scalar potentials and therefore we shall call them the gauge invariant potentials (g.i. potentials). 

It appears that one can completely rewrite the Lagrangian density ([l] - ^) using g.i. potentials according 
to (|l3|). We have then 



[Ft - VfY - (V X F) 



w 



q^F^^iVqf V + N{q^f + qt^)-2Mq^f, 



(14) 



where q^ = tpip and the remaining symbols have the meaning as before. 

Lagrangian density (|l]) - (^) is now automatically gauge invariant since it is defined by the gauge invariant 
quantities. Considering now F,/, q, n as new variables we obtain the system of equations 



Sf : (Ft - V/)j + V X V X F + 2q'^VF ^ 0, 
Sf:Af-V-Ft + 2q^ (M - Nf) = 0, 



Sn:Wn + Vn 



q^F^ + {Vqf - iV„ (q^f + g,^) + 2qHU = 0, 



VAq - Nqu -hWg~q (F^V - pN) 



FV-(V(Z)' Vg-N^ntqt 



Ng + Vn [Vn ■ Vq] - fq (2Af + qMg) = 



Rewriting the first two equations as before, we have the nonhomogeneous wave equations 



(15) 
(16) 

(17) 
(18) 



AF - Fu = 2gVF+V (V • F-ft) , 
A/ - ftt = -2g2 (M - iV/) + (V • F 
with the current and charge density defined as 

j = 2qVF, 

p - 2g2 (M - Nf) . 



(19) 
(20) 

(21) 
(22) 
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Now it is easy to check the continuity equations V -j+pt = 0. Indeed, by ( p^ ) and ( p^ ) 
V ■ j = V • {2qVF) = -VFtt + Aft = -2 [q^ (M - Nf)]^ = -p*. 

Two possible charge density components are given by ( p^ ) and it is only a question of a concrete physical 
problem which of them can be neglected, putting the relevant constant equal to zero. (The are no inevitable 
arguments to assume = 0; e.g. in paper there is considered the configuration when p ~ tp.) 

As in the case where on the vector and scalar potentials the Lorentz gauge could be imposed, now it can 
be also done formulating the Lorentz-like gauge imposed on g.i. potentials in form of condition 



V • F - /t = 0. 



(23) 



This requirement itself is also gauge invariant since it contains gauge invariant quantities, in contrast to 
the Lorentz gauge condition for A and / , which is not invariant with respect gauge (Q), unless AO — Ou — 0. 
It is obvious, since in the language of standard potentials equation ( p3[ ) reduces to V ■ A + tpt = AO — Ott- 

When the static solutions are considered f = ip and the system of equations (Esl) - (p^) reduces to 



AF - 2q^V'F = 0, 

A/ + 2g2 (M - Nf) = 0, 



VAq- 



q^Y^ + [Vqf - N„q^f + 2qHl„f = 0, 



(24) 
(25) 
(26) 



(27) 



(F^q^ - (V?)') V, - q^pN,'\ - q {F^V - pN) + K [Vti • Vg] 

where we assumed the static reduction of the Lorentz gauge for q.i. potentials V-F = 0. The current and charge 
density definitions ( [l5| ) and (^6|), respectively, remain unchanged. 

Defining the canonical momenta PF,P/,Pg,Pn associated with variables F,/, g, n in a standard form, we 

have 



PF = Ft-V/, p/ = 0, Pg = 2iVgt, p„ = 0. 
Therefore the Hamiltonian density TL — Ft ■ 'Pf ^ ItPq—J^ is 



(28) 



vl + {yxFY +V/-PF 



AN 



W 



q^F^ + iVq)' V + 2Mq^f ~ Nq^f 



(29) 



where the whole term [pq] /AN drops out if = and then also = 0. The standard Hamilton equations 
reproduce ( p8| ) and equations (|l^ - |l^), as should be. Substituting ( p8| ) into (p^), the energy density is 



F2 



t -(V/)' + (VxF)' 



W 



q^F^ + {Vqy 



V ■ 



(30) 



2Af + N 



iltf-f 



It is natural to ask about the conservation relations. Tedious, though simple calculations lead to the 
equation 



{Ftf + {rotF)' - (V/)' 



W 



q^F^ + {Vqy V + 2Mq\f + N {qtY - f 



V • (Ft X rotF+ {Ft - V/) U + 2VqtVq) 



(31) 



Till now, still the real functions W, V, M, N are completely arbitrary and some of them can vanish. Observe 
however, that the right hand side of ( |3l| ) i.e. flux depends only on V . Moreover, meanwhile we did not impose the 
Lorentz-like gauge (03) and also in particular cases the variable n can vanish simplifying all previous relations. 
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Admitting however the Lorentz-hke gauge (|2^) system of equations (|T^) - ( pO| ) becomes 



AF - Ftt = j = 2q^VF, 

A/ - ftt ^-P= -2q^ {M - Nf) , 



Nn {q^f + Qt^) + 2q^M,J = 0, 



VAq ~ Nqtt ^^Wq + q [F'V - fN) + ^ F'q' - {Vq) 



v„ 



(32) 
(33) 
(34) 

(35) 



q^f - {qtf Ng - Vn [WnVq] + N,,ntqt + fq (2M + qM^) , 



and in the next paragraph we shall discuss some particular situations due to the specific choice of W, V, M, and 
functions. 



3 Elementary reductions 

As a first example we consider the system of equations which follows from Lagrangian (|l|) - and from the 
relevant equations when the variable n does not appear at all, y = A^ = 1 and M vanishes. Moreover, assuming 
W {q) = q^ + /S^Z^/S, from (|3|) - ( |35| ) and (^ we obtain the system of equations: 



1 



{Ft - Vff - (V X F)' 



q'F' + {VqY +{q'f + qt') 



AF - Ftt = j = 2q^F 

A/ - ftt ^-P = 2q\f. 

"1 

2 



Aq-qtt = ^Wg + q{F^-f^ 
dt 



{Ftf + {rotFf - (V/)' 



W 



q^F' + {Vqy 



= V • (Ft X rotF+ (Ft - V/) ft + 2qtVq) 



(36) 

(37) 
(38) 

(39) 
(40) 



One can find elsewhere Its static version, but when the time dependence in included, it was discussed 
only in a few publications e.g. |^ fl^ ,. To our knowledge it has no broader application in superconductivity, 
probably because of two rather serious insufficiencies: the charge is proportional to the scalar potentials (see 
r.h.s. of ( ^sl ) and the second order derivatives with respect to time qtt in equation ( |39| ) appears. Rewriting this 
system in a standard language of Ajip^ip, we will have the second derivative iptt, which seems to be at least 
strange when the time dependent Ginzburg - Landau equation (or Gross - Pitaevskiii eqn.) are compared with 
the classical Schrodinger equation. 

In order to derive the system of equations with the first order derivative 4't in relevant GL equation, as e.g. 



in pap er ||4| one can consider a dual situation when V 
■"321) - (|3|)and (|l|) reduces then to 



M — I and A^ vanishes ( also n = 0). The system 



{Ft - Vff - (V X F)' 



AF - Ftt = J = 2q^F, 

A/ - ftt - -P = -2q^ 
1 
2 
'1 
2 



A? = ^VF, + qF^ + 2fq, 



dt 



{Ftf + {rotFf ~ {Vff 



W 



W 



g^F^ + {VqY 



q'F' + {VqY 



V • (Ft X rotF+ {Ft - V/) ft + 2qtVq) . 



2q'f, 



2q'f 



(41) 

(42) 
(43) 

(44) 
(45) 
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In contrast to the previous choice, now in language of g.i. potentials, equation (Q) does not contain the 
time derivative, but it is present when the system is written by means of A, if, ■0 variables. In fact, instead of 
(p4[), equation (^) reduces to form 

z9tV = ^(*V + A)'^ + iw^,+i/;^. (46) 

A comparison of both above discussed examples speaks in favor of the choice iV = 0, since it leads to the first 
order in time GL equation. 

In the next paragraph we consider a further argument in favor of that choice. 



4 Application to the superconductivity 

In a series of papers jl^, there was considered the Bardeen phenomenological extension of the static GL 
theory, mainly in order to include the energy of electric field and to demonstrate an important role of the scalar 
potential which is very often neglected. Unfortunately this theory is completely static, starting from the time 
independent free energy. Neglecting the unimportant details and in symbols (signature) adopted in present 
paper the free energy discussed there has the form 



(V X A)^ 



There is introduced a new variable n 

parameter. The function w includes the Gorter-Casimir corrections and w 



U + W+ — + A)^\'- + (47) 
4m 

the density of the total charge, while V' is reserved for the order 



lrp2 

2^ 



7- 



Essential for further consideration is that 7 = 7 (n) and thus w = w , n) . Moreover, the effective mass of 
the carriers also depends on density of charge i.e. m = m (n) . The function U introduced in ( ^7| ) described the 
effect of the screening on the Thomas - Fermi length and what is important U — U (n) and thus dU/dn = Cp, 
where C is constant. From the context of the paper it follows also that authors assume p ~ en. 

According to the authors, variations with respect A, Lp, ip and n variables lead to the system of equations 



VxVxA = Vxi7 = 

-Aip = p 

eifi = —Cp — w. 



2m 



Re 



1 1 

4 m 



(^-ip (iV + A)ip 



d (In m) 
d (Inn) 



(iV + A)V'| 



(48) 
(49) 

(50) 
(51) 



Accepting without discussion the physical aspect and phenomenological justification of the free energy 
according to expression ( ^7|) we shall show that one can propose an expression for the Lagrangian which includes 
time dependence and next to derive the self-consistent time dependent field equations which quite completely 
reduce to the system ( ^8|) - (js^) , when only static effects are considered. " Quite completely" means here that - 
up to the terms (V (n/m)) . The reasons that the terms proportional to V (n/m) can be neglected in equation 
( ^T|) was not specified. A proposition of a proper Lagrangian reduces in practice to a suitable choice of four 
functions W, V, M and N .in our considerations. 

Let assume that 



V 



v( 



2 IV'I 



TV 0, 



where still m = m (n) , w 



n/{4m), W = an^/2 + w, M: 
) ^IV'I^ and a is constant. Lagrangian (^) takes the form 



[At + v^y 



(VxA? -(an^/2 + w) 

J ^ ^ 4to 

an 



;iv + A)vr 



(52) 



(53) 
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It differs from the free energy ( |4^ ) practically only by time dependent terms. Consequently, the proper field 
equations (||) - (^) are 



(5a : AA - Att - — 

5^ : Aip - iptt - 



4m L V 

-an - (VA + ipt)f 



Sn a{n + ip) + Wn + \ ( |(iV + A)V'|^ 

4 \ni/ n 2 |-0| 



V (V • A + (^t) , 
(^tpdiip - ipdt-ip^ = 0, 



Tl - 



0. 



(54) 
(55) 
(56) 

(57) 



Assuming a Lorentz gauge VA + ift~Q and observing that (ji/m)^ — [I ^ d (Inm) /d (Inn)] /m, it is seen 
that static version of the system ( p^ - ( p7| ) is quite equivalent to system ( ^8| ) - (|5T|). The difference is only 

because of the term — -j ^-^Vt/) — iA IV'I^^ V {n/m) in the equation (|5^), which in (^ij) was probably omitted or 

simplified. 

We would underline that the time dependent system (|5j) - ( ^7| ) is self-consistent with all consequences 
which were mentioned in the previous paragraphs. The "prototype" of the GL equation ( ^ ) contains only the 
first order time derivative. 

Rewriting all equations in the language of the gauge invariant potentials F and /, the whole theory becomes 
much simpler, particularly when the Lorentz-like gauge is adopted. Substituting (^) to (^4|) - ( |l8| ) we obtain 
Lagrangian 



(Ft - Vf f - (V X F)' 



an- 
~2~ 



w 



Am 



and relevant field equations 



5f : AF - ¥tt a^F 



5/ : A/-^ + a7i= (V-F-/0,, 



J„ : an + Wn + ( t~ ) 
V 4to / n 

5„ : —Aq-2w„ qF^ 

m m 



-V (V ■ F - 



t)t ' 



q^F^ + iVqy +a/ = 0. 



Vn • Vg = 0, 



together with a conservation law (Poyinting like theorem) 



an/, 



(58) 

(59) 
(60) 
(61) 
(62) 



1 r 

2 



(Ftf + {rotFf - (V/) 
: V ■ f Ft X rotF+ (Ft 



n 

Am 



q'F^ + iVqy 



na 



(n + 2/) 



V/)/t 



2m 



(63) 



It is seen that the system ( p8| ) - (|6^) is much simpler than those expressed by the standard potentials, 
especially when the condition V • F — /j = is imposed. Moreover, the time derivatives do not appear at all in 
equations ( |6l| ) and (|6^), suggesting to use the static solutions as a first step when the perturbation technique 
will be applied for the solution of a dynamic problem. 

Repeating the calculations of this paragraph in case when iV ^ and it is e.g.N = f3n/2 , where /? is 
constant, the static versions of equations (|5J) and (|5^) do not change but the static version of equations ( ^5| ) 
and (^ ) take a form 



Aip = -n{a + (3ip); a (n + ip) + Wn + ] (—) | (iV + A) - ^</5^ 

4 Vm/ n 2 



= 0. 



(64) 



It is seen that the agreement of this system with the system (^8|) - (^ij) , which was proposed on a base of 
physical premises can be reached only when /? = 0, i.e. when in the starting Lagrangian (|) -(|) TV function 
vanishes. The proposed term, proportional to A/, will be the sole one which includes time dependence. 
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5 Conclusion 



We proposed a general and extended formula for the Lagrangian in spirit of the Ginzburg - Landau formalism, 
augmented by the time-dependent terms of two types. Moreover, a few arbitrary weight scalar functions 
introduced in this Lagrangian allow us to derive a family of the self-consistent time - dependent extended 
Ginzburg - Landau equations which coincide with the relevant Hamilton equations. Among them there are 
families with first and second time derivatives of the order parameter functions. The whole theory and equations 
can be simplified when it is rewritten by means of the gauge - invariant potentials. As an example of the 
application, the static version of GL equations reported in |0, are derived and verified as the reduction of 
a full dynamical theory. 

This reproducibility seems to indicate that for the correct dynamic theory, the scalar function TV discussed 
in the text should vanish. This means that then only the first order time derivative will be present in the time 
dependent version of GL formalism. 
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